Microscopic description of nuclei 
in the middle of the pf-shell 
by a shell model calculation 
with G-matrix interaction 

Hitoshi Nakada 

Department of Physics, School of Medicine, Juntendo University 
Hiraga-gakuendai 1-1, Inba-mura, Inba-gun, Chiba 270-16, Japan 

Takashi Sebe 

Department of Apphed Physics, College of Engineering, Hosei University 
Kajino-cho 3-7-2, Koganei, Tokyo 184, Japan 

and 

Takaharu Otsuka 

Department of Physics, Faculty of Science, University of Tokyo 
Hongo 7-3-1, Bunkyo-ku, Tokyo 113, Japan 



Abstract 



Energy levels and electromagnetic properties of 24 nuclides with AT = 28 ~ 30 are stud- 
ied in terms of a large-scale shell model calculation, which contains no newly adjusted 
parameters. The Kuo-Brown G-matrix interaction is shown to reproduce energy levels 
of 205 low-lying states of all these nuclei. We evaluate effective charges by incorporating 
the core-polarization effects caused by the coupling to GQR's. We then compute E2 mo- 
ments and transition probabilities. The Ml moments and transition rates are calculated 
by quoting the effective (^-factors of Towner, which are obtained by taking into account 
the meson-exchange and the core-polarization mechanisms. By this microscopic calcu- 
lation most of the E2 properties and the magnetic moments are reproduced. Although 
there are agreements and disagreements in the Ml transition rates, the general tendency 
is reproduced. The (e, e') and {PtP') excitation from the ground state to some low-lying 
2"*" states is also discussed. 



1 Introduction 



The nuclear shell model has been a basic tool to understand the low energy phenomena 
of nuclei in terms of the nucleonic degrees of freedom. Although the shell model has been 
successful in the p- and sd-shell regions[0, it is prohibitively difficult to apply the model 
to most of heavy nuclei, because of the large size of the model space. The effective inter- 
action between valence nucleons should originate in the interaction between free nucleons. 
There have been many investigations to derive an effective interaction from a microscopic 
viewpoint of this sort. The restricted model space which enables us to carry out the shell 
model calculation makes it necessary to renormalize the interaction. The renormalization 
appears to be successful in comparison to experiment, although the convergence of the 
renormalization procedure is still being investigated. Nevertheless, it is known that re- 
alistic interactions derived from free-nucleon interactions substantially explain the main 
features of low-lying states of sd-shell nuclei. An example is the interaction obtained in 
the earlier work of Kuo and Brown [0]. 

In the pf-shell, shell model calculations have been attempted for many years. In 
most earlier calculations, the excitation across the gap at Z, N = 28, as well as those at 
Z, N = 20 or 40, was ignored. Empirical effective interactions were often employed in 
those worksp, ^ Horie and Ogawa reported a shell model calculation for = 29 and 
30 nuclei|^, |^, which reproduces yrast levels systematically. A typical example is ^^Fe. 
We, however, cannot reproduce higher states in the framework of Horie and Ogawa's 
approach. For instance, the excitation energy of the 2j^ state of ^^Fe is overestimated by 
0.4MeV, and that of the 2'^ state by 0.5MeV. This discrepancy is serious in discussing 
the mixed-symmetry collective states|H], which have been attracting a certain amount of 



interest §, |, |Tg, [n|, |T2[ . 

In contrast to the phenomenological interactions, Kuo and Brown applied their G- 
matrix method to the pf-shell region, and derived a realistic interaction on the top of the 
^°Ca core[|l^]. The Kuo-Brown interaction has been applied to shell model calculations in 



the full pf-shell, for 40 < A < 44 nuclei |T^. Though qualitative agreement with the exper- 
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imental data is found in those calculations, there still remain recognizable discrepancies. 
A few ways of modification of the interaction have been proposed [p!5| , |T6i [T7| . However, 
we should take notice that excitations from the sd-shell may be significantly present in 
the mass region of 40 < A < 44. It is worthwhile to test the Kuo-Brown interaction in 
a sufficiently large model space around the middle of the pf-shell, where low-lying states 
will minimally be affected by other shells. 

Oberlechner and Richert carried out a shell model calculation around ^^Ni|jl8|], by 
using another realistic interaction. They, however, could not succeed in reproducing the 
experimental energy levels so well, presumably because of the truncation of the model 
space. 

There have been several attempts to take into account the one-particle excitation 



from the Of 7/2 to the IP3/2 orbit [|T9l. There have been shell model calculations which 
included the one-particle excitation from the 0f7/2 to any of the 0f5/2, IP3/2 and lpi/2 
orbit s [^ . However, the situation has not necessarily improved by this extension of the 
configurations; the discrepancies with Ex{2f) and Ex{2f) of ^^Fe, for instance, remain 
similar to the results in Ref.|0]. This probably happens because the excitation of two 
nucleons from the Of 7/2 to the higher orbits is crucial and can indeed be driven by the 
deformation and/or the pairing correlation. 

In this paper we show the results of the large-scale shell model calculation for 20 < Z < 
28, 28 < < 30 nuclei, including the excitation of up to two particles (four for ^^Ni) from 
the Of 7/2 to the other pf-shell orbits, which seems to be of particular importance in this 
region, as mentioned just above. It turns out that, if a sufficiently large configuration space 
is chosen, the Kuo-Brown interaction appears to be quite an excellent effective interaction. 
The model space and the effective hamiltonian are discussed in more detail in Section |^. 
In Section ^ calculated energy spectra are shown, in comparison with the observed ones. 
In Sections ^ and ^, the E2 and Ml properties of those nuclei are analyzed in terms of the 
operators obtained microscopically. Those results will tell us about the reliability of the 
shell model wavefunctions, and also how well we can understand the nuclear properties 
from the nucleonic degrees of freedom. The (e, e') form factors concerning the quadrupole 
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collectivity are discussed on the same footing as the E2 transitions, in Section |^. The 
DWBA results of (p, p') differential cross sections are briefly surveyed in Section 0. They 
also supply a justification of the shell model wavefunctions relevant to the low energy 
quadrupole collective motion. Section ^ is a summary of this article. 



2 Model space and effective hamiltonian 

We first show the configuration space of our shell model calculation. Assuming ^°Ca to 
be a doubly magic inert core, we consider configurations such as 

(0fr/2)"^-^0f5/2lP3/2lPl/2)"^+', (1) 

where rii and 77-2 are defined so that k = should give the lowest configuration in which 
the excitation across the gap between 0fi/2 and the other three orbits is absent. For 
instance, rii = 14 and n2 = 2 for ^^Fe, since this nucleus has 6 protons and 10 neutrons 
in the pf-shell. For the nuclei handled in this paper, rii = {Z — 20) + 8, ^2 = — 28. The 
present model space contains all configurations of /c = 0, 1 and 2 for each nucleus. We 
would like to emphasize that such a large configuration space has never been used except 
for Refs.p, [l^, ^ The k = 2 configurations evidently play an essential role in the 
description of the relaxation of the magicity of Z = A = 28 due to the pairing correlation 
and the deformation. 

Concerning the effective hamiltonian, we adopt the Kuo-Brown hamiltonian on top of 
the ^°Ca core [13]. In this hamiltonian, the single-particle energies have been determined 
from experiment; eof,/^ = O.OMeV, cof^/^ = 6.5MeV, ei^^^^ = 2.1MeV and eip^^^ = 3.9MeV. 
The two-body interaction is derived from the G-matrix based on the Hamada- Johnston 
potential, assuming the single-particle wavefunction in the harmonic oscillator approx- 
imation. We include the 3p-lh correction, which is evaluated by the perturbation. It 
is noted that there are only few reports of shell model calculations in the pf-shell using 
such a microscopic interaction without any phenomenological modification. Around the 
middle of the pf-shell, both the influence of the sd-shell and that of the sdg-shell are ex- 
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pected to be minimal. The Kuo-Brown interaction seems to be valid in this region, if the 
configuration space is sufficiently large. Indeed, the success of the Kuo-Brown interaction 
in the k < 2 configuration space was reported for ^^Sc in Ref. |T^. We already showed the 
validity also in ^^FeQ], and the extensive application is carried out in this study. 

The isospin is automatically conserved, and in the following we will deal with states 
with the lowest isospin for each nucleus. 

In the numerical calculation, we use an M-scheme shell model code newly developed. 
The Lanczos method is applied for diagonalizing the shell- model hamiltonianj^. The 
present model space is one of the largest one available at present. In the M-scheme, the 
M = space has the largest dimension in even-mass nuclei, while the M = ±i space in 
odd-mass nuclei. This dimension is displayed for each nucleus in Table as well as the 
largest dimension in the JT-scheme. 

We add some remarks on the effective interaction: the realistic interaction has long 
been discussed in the sd-shell. The bare G-matrix is insufficient to reproduce the ob- 
served energy spectra, and some renormalization is required. Though the 3p-lh correc- 
tion improves the energy spectra significantly^, the convergence of the renormalization 
procedure has been a serious problem. Despite much effort and a certain progress ||24]|, 



this problem has not been fully solved yet, and it is beyond the scope of this article. We 
only comment that the 3p-lh correction gives reasonable spectra and that the low-lying 
energy spectra resulting from an interaction including the higher-order corrections do not 



look quite different from those obtained by the 3p-lh correction p5|. 

In the 20 < Z < 28, 20 < < 28 region, there are many intruder states, which are 
not described within the shell model for the pf-shell. These intruder states will primarily 
consist of the excitation from the sd-shell orbits to Of 7/2- Whereas one hole in the sd-shell 
leads to negative-parity states, two-hole configuration could be low in energy owing to 
collectivity. The proton-neutron interaction should play a dominant role in introducing 
this collectivity. The coexistence of the intruder configuration significantly infiuences the 
spectrum, and cannot be treated within the renormalization scheme of Kuo-Brown's 3p-lh 
correction, which is based on the perturbation theory. The modification of the Kuo-Brown 
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Table 1: 

The largest dimensions in the shell model calculation with k <2 configurations, both in the 
M-scheme and the JT-scheme. The fourth and fifth columns indicate the spin-parity and 
isospin for which the dimension in the JT-scheme is largest. 



N 


nucl. 


M-scheme 


J 


T 


JT-scheme 




49Sc 


1,561 


7 9 
2' 2 


7 
2 


210 




50xi 


5,530 


4 


3 


656 




Sly 


11,676 


9 
2 


5 
2 


1,254 


28 


52 Cr 


16, 544 


5 


2 


1,615 




53 Mn 


15,883 


9 
2 


3 
2 


1,456 




54Fc 


10,620 


4 


1 


891 




55Co 


4,717 


9 
2 


1 
2 


347 




56Ni 


1,353 


3 


1 


103 




50Sc 


5,550 


4 


4 


740 




5iTi 


20, 425 


9 
2 


7 
2 


2,423 




52y 


46,034 


5 


3 


4, 963 


29 


53 Cr 


68, 355 


9 
2 


5 
2 


6,929 




54 Mn 


70, 740 


5 


2 


6,751 




55 Fe 


50,425 


9 
2 


3 
2 


4,611 




56Co 


24, 844 


4 


1 


2,126 




57Ni 


7, 890 


7 
2 


1 

2 


619 




51Sc 


13,411 


9 
2 


9 
2 


1,774 




52Ti 


52, 624 


4 


4 


6,207 




53v 


123,519 


9 
2 


7 
2 


13, 574 


30 


54Cr 


195,334 


5 


3 


20, 074 




55 Mn 


211,833 


9 
2 


5 
2 


20,992 




56 Fe 


162,358 


5 


2 


15,457 




57Co 


84, 978 


9 
2 


3 
2 


8,029 




58Ni 


29, 792 


4 


1 


2,780 
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interaction discussed in Refs. |T5|, |T6[ yields a better agreement in the 20 < Z < 28, 20 < 



< 28 region, and this may be connected with the influence of the intruder configuration. 
The influence of the intruder configuration becomes weaker in > 28 nuclei, due to the 
suppression of the neutron excitation from the sd-shell to Of 7/2 by the Pauli blocking. 
This will also reduce the proton excitation from the ^°Ca core through the proton-neutron 
correlation. On the other hand, as we approach Z = N = 28, the excitation from 0f7/2 
to (0f5/2lp3/2lPi/2) "will grow. In such cases, the extended configuration space currently 
considered becomes crucial, as can be seen in the following sections. 

In Ref. |T^, Kuo and Brown tabulated the interaction for the model space consisting 
of the four pf-shell orbits and Ogg/2- In this study we omit the contribution containing the 
0g9/2 orbit, because the influence of 0g9/2 is not expected to be significant for the nuclei 



under consideration. We will return to this problem of the model space in Sections ^ 
and |[ 

We shall use a mass-number independent interaction; the harmonic oscillator bases 



have been employed in calculating the G-matrix, with hcu = lO.SMeV in Ref.[|13|. The 
oscillator length assumed usually as 5 oc A^^^ varies by only a few percent among the nuclei 
under consideration, namely A = 49 ~ 58. The above value hu = 10.5MeV adopted in 
Ref. HTSf corresponds to A ~ 60 if we take hu = 41.2^4"^/'^. Thus, the G-matrix calculated 
in Ref. [O] appears to be appropriate for the mass region considered in this study. 



3 Energy levels 

In this Section we systematically show results of the present calculation for energy levels 
of 20 < Z < 28, 28 < < 30 nuclei, in comparison with experimental data. The data are 
quoted from the compilation in Refs.||26|-|35[|. We concentrate on the excitation energies 



and do not discuss the binding energies. 
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3.1 Ground-state configuration and convergence by truncation 

Before going to the energy levels, we shall see the ground-state configurations classified 
in terms of k in Eq.(|l|). They are shown in Table |^. 

Though one might be interested in the Ca isotopes, they are not suitable for the 
present study. Probably, in the Ca isotopes, we need to take into account the excitation 
from the sd-shell. 

Within the present calculation, it turns out that the ground state of any nucleus is 
dominated by the k = configuration, as is expected. The ground state of ^^Ni is outside 
the (0f7/2)^^ closed-core by 37%. The k = probabilities are around 60% except for Sc 
isotopes, while about 70% in Sc isotopes. Thus, the ground state is dominated by the 
k = configuration in all cases studied. The importance of the k = 2 configuration is 
obvious in Table Q The probability of the k = 2 configuration is even larger than that of 
the k = 1 configuration in many of the nuclei under study. 

Table ^ shows that, for = 28, 29 and 30, the k = 2 probability goes up as Z 
increases. This observation gives us information whether the wavefunction is converging 
under the k < 2 truncation. This convergence will be referred to as /^-convergence. As 
Z is larger, the proton Fermi energy should be higher. Therefore the valence proton is 
more susceptible to be excited from Of 7/2 to the upper orbits. This proton excitation will 
induce the excitation of neutrons from Of 7/2 through the proton-neutron correlation. In 
other words, the energy gap above the 0f7/2 orbit becomes less important. Consequently, 
as we approach Z = 28, the A;- convergence would become relatively worse. The sd-shell 
configurations, on the other hand, could have a stronger infiuence for Z close to 20, since 
the Pauli blocking suppresses to a lesser extent the proton excitation from the sd-shell to 
Of 7/2. These problems seem to be much less serious in the region Z ~ 24. 

The k = 2 probability is slightly smaller in = 29 and 30 than in A^ = 28 isotones. 
This happens probably because, in A^ = 29 and 30 isotones, the partial occupancy of the 
neutron (0f5/2lp3/2lPi/2) orbits suppresses the excitation from 0f7/2. Taking into account 
the diminishing excitations from the sd-shell, the /^-convergence is expected to be quite 
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Table 2: 

Ground-state configuration in the shell model calculation in terms of the probabilities (%) of 

each k configuration. 



N 


nucl. 


state 


k = 


k = l 


k = 2 


28 


49 SC 
50Ti 

Sly 
52Cr 
53 Mn 

55 Co 
56Ni 


of 

[|]r 


68.2 
61.1 
60.5 
60.8 
61.3 
62.2 
63.0 
63.4 


10.3 
14.9 
14.7 
12.3 
10.8 
7.0 
4.5 
0.0 


21.5 
24.0 
24.8 
26.9 
27.9 
30.8 
32.5 
36.6 


29 


50Sc 
51xi 
52v 

53 Cr 

54 Mn 
55Fe 

56 Co 

57Ni 


[fir 
ir 
ir 
[fir 


68.8 
60.7 
56.6 
57.3 
55.5 
58.4 
59.4 
60.9 


14.6 
19.2 

22.9 
20.2 

21.6 
15.1 

12.1 

6.3 


16.6 
20.1 
20.5 
22.5 
22.9 
26.5 
28.5 
32.8 


30 


51SC 

52Ti 
53y 

54 Cr 
55Mn 

57Co 
58Ni 


[|]r 
[fir 
[|]r 


70.2 
62.7 
59.8 
56.5 
53.2 
55.7 
56.8 
57.7 


10.8 
14.8 

19.3 
20.6 

25.1 
17.7 

15.4 
9.1 


19.0 
22.5 
20.8 
22.9 
21.7 
26.6 
27.9 
33.1 
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good for = 29 and 30 nuclei. 

According to the above considerations, we anticipate that, among the nuclei under 
discussion, the /c-convergence will be worst for ^^Ni. On the other hand, it is possible 
to diagonalize the shell model hamiltonian for ^^Ni in the space containing all the k < 4 
configurations, which leads to the M-scheme dimension of 497, 805. By this calculation, 
the ground-state energy of ^^Ni is lowered by 3.3MeV, and the k = probability reduces 
to 28%. We shall return to this k < 4 calculation in Subsection |3^ . 

While there seems to be a slight difference in the fc- convergence of the wavefunctions, 
the energy levels are excellently reproduced, as is discussed in detail in the subsequent 



subsections. We will discuss the convergence also in Sections H and KB 



3.2 ^^Ni, 49Sc and ^^Co 

In the A; = configuration space, we have a single state with = 0+ in ^^Ni, and one 
with = \ in ^^Sc and ^^Co. These spin-parities of the ground states are unchanged 
in the present calculation, in agreement with experiments. 

The first excited levels become too high in ^^Ni and ^^Co in the k < 2 calculation. In 
order to describe these states, one probably has to consider a softer ^^Ni core than the 
present one. Higher k configurations will be necessary in order to reproduce those levels. 

Wong and Davies carried out a shell model calculation for ^^Ni by including the k < A 
configuration in part 1 22]. They reported that the Kuo-Brown interaction successfully 



reproduces the energy spectrum of ^^Ni, although their single particle energies are some- 
what different from the present ones. There the (0f7/2)^^ probability of the ground- 
state wavefunction was 16%. The small (0f7/2)^^ probability has brought about several 
discussions |15|, |T8|. In the k < A calculation for ^^Ni with the present hamiltonian, the 
ground state seems to have a similar structure to the result of Wong and Davies, though 
the k = probability is larger in the present result. In Table |^, the excitation energy of 
the 2f level obtained by the k < 4 calculation, as well as by the k < 2 calculation, are 
shown in comparison with the datum. The k < 2 calculation yields too high an excitation 
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Table 3: 

Comparison of the k < 2 and k < A calculation with the experimental data taken from 



Refs-lgg, for Ex{2"[) (MeV) and ■B(E2;2]^^ 0^) values (e^fm^) in ^^Ni 



quantity 


Cal. {k < 2) 


Cal. (A; < 4) 


Exp. 


Ex{2t) 


[MeV] 


5.792 


2.724 


2.701 


i?(E2;2+^0+) 


[e^fm^] 


80.4 


80.6 


74 ± 11 



Figure 1: 

Experimental and calculated energy levels of ^^Ni. The experimental data are taken from 



Ref.|34]. The calculated energy levels are obtained from the present k < 4 shell model 

calculation. 

energy, suggesting that the k < 2 model space is too small for this state. On the other 
hand, the observed energy is reproduced quite well by the k < A calculation. This agree- 
ment encourages us in applying the present G-matrix interaction to the nuclei around 
^^Ni. The k < 4 calculation yields a O"*" level at Ex = 2.30MeV, which is strongly domi- 
nated by the k = 4 configuration. An unobserved 2+ level also appears at Ex = 2.96MeV 
in this calculation. It is interesting to know whether or not these levels exist. The low- 
lying levels of ^^Ni obtained by the present k < 4 calculation are summarized in Fig.[l|, 
in comparison with the data|^. We note that the indispensable roles of /c = 3 and 4 
configurations in ^^Ni are due to Z = = 28 and also due to that neither excited 0^ nor 
2"*" is allowed with A; = 0. 

In ^^Sc, no negative-parity states have been observed below Ex = 3MeV other than 
the ground state. Nevertheless there emerge about ten negative-parity states in 2.4 < 
Ex < 3MeV in the calculation. Since ^^Sc is far from /3-stable line, those states may not 
have been observed yet. Because Z ~ 20, the influence of the excitation from the sd-shell 
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Figure 2: 

Energy levels of ^''Ni. The experimental data are taken from Ref.[^^. The calculated energy 
levels are obtained from the present k < 2 shell model calculation. 

may be sizable in comparison with the larger Z case. 



3.3 5^Ni 



Suppose that Ni is taken to be an inert core, there is only one valence neutron in the 
^^Ni nucleus. Therefore we have only three states with J"^ = | , | , | in the A; = 
configuration. 

In Fig.|^, the calculated and measured energy levels are compared. The agreement is 
sufficiently good for the lowest three states, which predominantly consist of the k = 
configuration. On the other hand, the energy gap between the lowest three states and 
the other states is larger in the calculation than in the experiment, similarly to ^^Ni and 
^^Co. A larger model space will be necessary for these higher states. 

In the k = space with the same hamiltonian, we have too low an energy for the 
state {Ex = 0.28MeV). The k = 1 and 2 configurations improves the energy by raising 
the relative to the and states. 

Note that no = I"*" state is observed as low as the , [|]r £^iid [■^ji levels. A 
candidate of [Ij^*" is reported experimentally around Ex = 3MeV at lowest p4| . This fact 



suggests that we do not have to include the Ogg/2 orbit explicitly, as far as we restrict 
ourselves to the states below about 3MeV. 

Recently phenomenological interactions in the pf-shell, named FPD6 and FPMI3, 
have been proposed|0. These interactions are adjusted for the full pf-shell calculation 
in 40 < A < 44 nuclei and Ca isotopes. It is worthwhile testing these interactions in the 
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Figure 3: 

Energy levels of ^''Ti. The data are taken from Ref . [p7| . 



Figure 4: 

Energy levels of ^^Cr. The data are taken from Ref.|2S]. 



present mass region. 

In the k < 2 space, neither FPD6 nor FPMI3 reproduces the ^^Ni spectrum; FPD6 
makes the lower than the [^]{ by O.SMeV, while in the case of FPMI3, both the 
and the become lower than the by 2.5 and 1.4MeV, respectively. Note that 
the probabilities of the k = configuration are larger than in the case of the Kuo-Brown 
interaction for those three lowest-lying states. Thus, those FPD6 and FPMI3 results are 
less sensitive to higher k configurations. Both FPD6 and FPMI3 are work better for 
lighter pf-shell nuclei | 



3.4 Z = even, N = 28 nuclei 

The energy levels in ^''Ti, ^^Cr and ^^Fe are displayed in Figs.|^, ^ and ^. We find that the 
low-lying states correspond well to the observed levels, in the energy range of Ex ^ 4MeV, 
and with deviations of 6E ^ 0.3MeV. 

In 5°Ti, the 4^ and 3^" states appear around 3.5MeV in our calculation, whereas they 
are not seen in experiments. The occurrence of such additional levels resembles the case 
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Figure 5: 

Energy levels of ^^Fe. The data are taken from Ref . []3l| . 

of "^^Sc. The calculated 0^ state corresponds well in energy to the unestablished 0+ state. 

In ^^Fe, the number of observed states is larger than that of the calculated ones, 
probably due to the model-space restriction. We note that the observed and 22 states, 
which are missing in the present calculation, might belong to a deformed rotational band, 
since the energy interval between the two states is much smaller than that between Of 
and 2t. 

Within the k = configuration space, we have only four states with = 0+, 2+, 4+, 6+ 
in ^°Ti and ^^Fe. These states consist in the identical energy spectra between ^^Fe and 
^°Ti, as a consequence of the particle- hole symmetry. It is noted here that, within the 
k = configuration, Ex{2f) = 1.02MeV is obtained by the present hamiltonian for 
^°Ti and ^^Fe. The experimental energies are 1.55MeV for ^°Ti and 1.40MeV for ^"^Fe, 
respectively. The present calculation within the k < 2 space improves the 2f levels to a 
great extent, somehow overshooting. Namely, due to the additional pairing correlations 
associated with the k = 1 and 2 configurations, the 2f energies are raised relative to the 
ground-state energies. 

There are two 4+ states around 2.5MeV in ^^Cr. This approximate degeneracy is 
reproduced in the present calculation, where both of the two states predominantly consist 
of the k = configuration. The correspondence of these states between calculation and 
experiment is tentatively made according to the electromagnetic properties, as will be 
stated in Sections H and ^ The observed and 22 states might be deformed states, as 
in the case of ^^Fe. 
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Figure 6: 

Energy levels of ^^V. The data are taken from Ref . p8| . 



Figure 7: 

Energy levels of ^^Mn. The data are taken from Ref.|3C]. 



3.5 Z = odd, N = 28 nuclei 

As shown in Figs.^ and ^ the present calculation reproduces the observed levels in Z = 
odd, N = 28 nuclei with the energy range of Ex < 3MeV, which is smaller than in 
the Z = even case. This is a reasonable consequence, because ground-state energies in 
even-even nuclei are systematically lower than those in neighboring odd nuclei, owing to 
the pairing correlation. Indeed, in the present case, this is reflected in the higher level 
densities in odd nuclei around the ground states. 

By comparing Figs.§ and ^, a notable similarity is found for the lowest five levels 
between ^^V and ^^Mn, analogously to the similarity between ^°Ti and ^^Fe. Within the 
space consisting only of the A; = configuration, this can be understood as a result of the 
particle- hole symmetry. The results in Figs.^ and |^ include the A; = 1 and 2 configurations, 
which break the simple picture of the particle-hole symmetry. Nevertheless, for the lowest 
five levels, the properties of this particle-hole symmetry remain to a certain extent in the 
present calculation. Such properties are seen experimentally also. On the other hand, the 
similarity disappears in higher states with Ex > 2MeV. It is noticed that our calculation 
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Figure 8: 

Energy levels of ^^Ti. The data are taken from Ref . [p8| . 



Figure 9: 

Energy levels of ^^Cr. The data are taken from Ref.|3C] 



reproduces several levels higher than 2MeV. 



3.6 Z = even^ N = 29 nuclei 

The energy levels of ^^Ti, ^^Cr and ^^Fe are shown in Figs.|^, |^ and |T0|, respectively. The 
good agreement between the calculated and observed ones is confirmed for Ex ^ 3MeV, 
as in the case of the Z = odd nuclei discussed in Subsection |375| . 

In ^^Ti there is an unobserved | level at 2.9MeV (labeled as [|]^ in Fig.^) in the 
present calculation. Concerning this level, we find a similar pattern in the energy spectra 



Figure 10: 

Energy levels of ^^Fe. The data are taken from Ref.[32| 
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Figure 11: 

Energy levels of ^''Sc. The data are taken from Ref.[p7| 



of ^^Cr and ^^Fe, in which the [|]2" is observed just above the [|]|f. When the good 
agreement around this level and the systematics over these three nuclei are considered, 
the occurrence of the [|]2 state is very likely. It is of interest whether this level is observed 
in the future experiment or not. 

We now compare the present results with those by Horie and Ogawa]^, which have 
been the most successful shell model calculation in this region. It should be recalled 
that only the k = configuration is included in the Horie-Ogawa calculation. Low-lying 
levels in = 29 nuclei are fitted in determining the proton-neutron interaction in the 
Horie-Ogawa calculation. Up to Ex = 2MeV, the two results are very similar except for 
the [|]2 state, in any of ^^Ti, ^'^Cr and ^^Fe (See Figs. 3-5 of Ref. 0). The Horie-Ogawa 
calculation does not produce [|]^ levels in ^^Cr and ^^Fe, which have been observed close 
to in energy. This is consistent with our observation that the k = 1 configuration 
dominates the [|]^ state of ^^Cr and ^^Fe. The correspondence between the calculated 

2 

basis of electromagnetic properties. 
3.7 Z = odd, N = 29 nuclei 



and experimental ^ levels which are nearly degenerate is made in Fig.pl and [TO on the 



In FigsjTl|^I^, we see a good agreement between the calculated energy levels and exper- 
iments for the Z = odd, N = 29 nuclei, as far as the energy range of Ex ^ 2MeV is 
concerned. This smaller energy range of agreement than that for even-even or odd nuclei 
seems to be associated with the pairing correlation, as has been stated in Subsection [STS] . 
Suppose that we restrict the proton and neutron valence orbits only to 0f7/2 and IP3/2 
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Figure 12: 

Energy levels of ^■^V. The data are taken from Ref.|29|. 



Figure 13: 

Energy levels of ^^Mn. The data are taken from Ref.|31] 



Figure 14: 

Energy levels of ^^Co. The data are taken from Ref.|33| 
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respectively, only four states are possible in ^°Sc; = 2+, 3+, 4+ and 5"*". Indeed, the 
lowest four states consist predominantly of that configuration in the present calculation. 
The observed 0.257MeV level is expected to have = 2+. There are several states 
without the corresponding observed levels in ^°Sc. This problem could be accounted 
similarly as ^^Sc is explained. 

Among all nuclei that we discuss in this paper, the ^^V nucleus is the only one for 
which the ground-state spin is not reproduced. This may not be a serious drawback, since 
the lowest three states are very close in energy and such approximate degeneracy itself 
is reproduced very well. Based on the present calculation, the observed 0.017MeV state 
appears most likely to be a 2"*" state. The tentative assignment of to the 0.023MeV 
state seems plausible. 

In the Horie-Ogawa results (See Fig. 7 of Ref.0), the 4^ level has energy very close to 
the lowest three levels (5^^, 3f and 2'1), in contrast to experiment. The present calculation, 
on the other hand, produces a level nearly degenerate with 1^ rather than the lowest 
three levels, in agreement with the data. The level, which seems to correspond to 
the observed state at Ex = 0.881MeV, is dominated by the A; = 1 configuration. The 
occurrence of this level is one of the advantages of the present large-scale calculation, 
while this k = 1 dominance explains why the Horie-Ogawa interaction is not capable of 
describing this state. 

It is demonstrated that the lowest five states of ^'^Mn are excellently reproduced. We 
confirm that the present result is better than the result of the Horie-Ogawa calculation 
even for these lowest states. There are several k = 1 dominant states among the displayed 
levels; for instance, 6i, li, 22 and Tf. The low energy of the cannot be described by 
the Horie-Ogawa A; = calculation (See Fig. 8 of Ref. [^]), which yields Ex{6f) ~ 2MeV. 
Despite a certain underestimation (~ 0.4MeV), the reproduction of the low Ex{6f) is a 
prominent advantage of the present approach. 

In ^^Co, the energy of the level shows a certain improvement from the Horie-Ogawa 
result (See Fig.9 of Ref. 0). There is an unobserved 2+ state around Ex = 1.3MeV in 
the Horie-Ogawa calculation. Such a state is not seen in the present result. The 1^ state 
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Figure 15: 

Energy levels of ^^Ti. The data are taken from Ref.[p9| 



Figure 16: 

Energy levels of ^^Cr. The data are taken from Ref.|31]. 



has the k = 1 dominance in the present calculation, consistently with the failure of the 
description of this state in the Horie-Ogawa calculation. 

We cannot reproduce the observed low energy of Oi in ^^Co. We should remember, 
however, that there is no 0"*" state in the A; = configuration. Our calculated 0^ level 
at Ex ^ 2MeV, whose appearance is encouraging to us, may come down when higher 
configurations are included. 



Figure 17: 

Energy levels of ^^Fe. The data are taken from Ref.[33|. 
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Figure 18: 

Energy levels of ^^Ni. The data are taken from Ref.[]35| 



3.8 Z = even, = 30 nuclei 



As is already reported in Ref. p] and is demonstrated in Fig.^, the spectrum of ^^Fe 
seen in experiments is excellently reproduced for Ex <, 4MeV, apart from the 3~ state 
of 3.07MeV, which is obviously outside the present configuration space. Discrepancies in 
the excitation energies are less than 0.2MeV. The calculated excitation energies of the 
yrast states also agree with the experimental ones up to = 8^; 4.75MeV (4.70MeV in 
experiment) for 7^, 5.32MeV (5.26MeV) for 8^ 

The energy levels of ^^Ti and ^^Cr are displayed in Figs.^ and |16|, respectively. There 
are a few calculated but unobserved levels in ^^Ti. It is not clear, as in the case of ^^Sc, 
whether those levels emerge because the effect of the sd-shell configuration is too small 
in the present calculation, or whether they have not been observed. 

The comparison with the data of ^^Ni is presented in Fig.^. The calculated 0^ state 



is much lower than the experimental one. However, the calculated is highly domi- 
nated by the k = 2 configuration. It is possible that this remarkable difference of the 
wavefunction from the ground state makes it difficult to observe the state. If we consider 
the calculated state to correspond to the observed second 0"^, the agreement becomes 
satisfactory. We then predict a 0^ state around 2.2MeV. The 2f state is also dominated 
by the k = 2 configuration. The i?(E2; 2j^ 0^) value is fairly large (200[e^fm^]) in 
the present calculation, suggesting a quasi-band structure. Although, in Fig.|l8|, a corre- 
spondence of the 23" state to the observed 2^ state is tentatively indicated based on their 
energies, this correspondence can be reconsidered. The calculated 23" state might also 
remain unobserved. The comparison between calculation and experiment, however, is not 
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Figure 19: 

Energy levels of ^^Sc. The data are taken from Ref . [p8| . 



Figure 20: 

Energy levels of ^^V. The data are taken from Ref. pO|| . 
disturbed by this open question so much. 

3.9 Z = odd, iV = 30 nuclei 

The calculated and observed energy levels are exhibited in Figs.|T9|-^, for Z = odd, 
N = 30 nuclei. There are plenty of low- lying states whose spins and/or parities are not 
assigned, in these nuclei. We can give a theoretical suggestion. 

In ^^Sc, the spin-parity of the fourth lowest state is | in the calculation, while | < 
J < I is suggested in the experiment, based on the log ft value of the /3-decay from 



Figure 21: 

Energy levels of ^^Mn. The data are taken from Ref.|32] 
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Figure 22: 

Energy levels of ^^Co. The data are taken from Ref. 

^^Ca[^. It is desired to specify the spin-parity of this state. According to the present 
calculation, the observed state with Ex = 1.394MeV should have J"^ = | , which is 
consistent with the experimental indications. The spin-parity of the 1.715MeV state is 
considered to be | . 

In ^^V, the calculation strongly suggests the = | assignment for the Ex = 
1.266MeV state, |" for the 1.550MeV state and |" for the 1.653MeV state. None of 



these suggestions contradict the experiments |30| . 

In ^^Mn, the data imply an extraordinary situation that three states are present around 
Ex = 1.29MeV with surprising degeneracy {AE = O.OOSMeV). Only an y state appears 
in the calculation. 

The existence of the | state at 1.757MeV in ^^Co is not explained within the present 
model space. This state might be an intruder state dominated by A; > 2 configuration. 

As a whole, the energy levels of the Z = odd, N = 30 nuclei are pertinently described 
by the present calculation, up to Ex ~ 2.5MeV. 



3.10 Summary of energy levels 

We summarize this section. The energy levels of the 20 < Z < 28, 28 < < 30 nuclei are 
well reproduced by the present work, where the large-scale shell model calculation in the 
k < 2 space is carried out by adopting the Kuo-Brown interaction. The upper bounds of 
the agreement with the experiments are Ex ~ 4MeV for even-even, 2.5MeV for odd-mass 
and 2MeV for odd-odd nuclei. The discrepancy is 5E < O.SMeV. 

Most of the low-lying states are, as are expected, dominated by the k = configuration. 
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In odd-odd nuclei, a certain improvement from the Horie-Ogawa result, which contains 
only the k = configuration, is achieved even for the lowest-lying k = dominant states. 
Some low-lying states are found to have larger probabilities of the k = 1 configuration 
than of the A; = configuration; 3+ and 4^ of ^^Cr, , IyI^ and of ^^Mn, one of 
the nearly degenerate |~ of ^^Ti, ^^Cr and ^^Fe, 5^, 3^ and 4^ of ^^V, 6^, 1^, 2^ and 7| 
of 54Mn, 1+ of 56Co, 4+ and 2+ of ^^Ti, 2+ and 3^" of ^^Fe, 4+ of ^^Ni, [|]r, [i]^ and 
of ^''Co. These levels show good agreement to experiments. We can expect that the 
wavefunctions of the above k = 1 dominant states are as good as those of other low-lying 
states, since the k = 2 probabilities of such states are similar to those of the lowest-lying 
k = dominant states. It is remarked that these k = 1 dominant states have been beyond 
the description by the Horie-Ogawa calculation. 



By comparing several effective interactions derived from the G-matrices[0, |T^, it 
is inferred that the G-matrix interaction seems to involve an uncertainty of a few hundred 
keV. Since this error in the hamiltonian is transferred to the results of the diagonalization, 
the present discrepancy, 6E < 0.3MeV, appears to be quite reasonable. Further investiga- 
tions will be required in order to discern whether the success of the present calculation is 
accidental or not. A test of the wavefunctions through the observation of electromagnetic 
properties is particularly important. 

The occurrence of the O"*" level around Ex ~ 2.2MeV in ^^Ni appears to be a challenging 
problem. It is strongly desired to search for such a 0+ level experimentally. A search for 
a O"*" level around Ex ~ 2.3MeV in ^^Ni, which emerges in the k < 4 calculation, is also 
desired. 



4 E2 properties 

By using the shell model wavefunctions we can investigate the electromagnetic properties 
of the nuclear states. In this and the next sections we discuss the E2 and Ml properties 
of the = 28 ~ 30 nuclei. The purpose of this study is, as well as to test the shell model 
wavefunctions, to see to what extent we understand the electromagnetic properties of the 
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nuclei in this region from the nucleonic degrees of freedom without introducing adjustable 
parameters. We employ single-particle parameter sets of the transitions derived from 
microscopic theories, which seem to be the most plausible ones presently available. We 
do not attempt to improve those theories nor to adjust the parameters in this work. The 
experimental data are taken from Refs.p6|-lB5|. 



4.1 Single-particle parameters 

The E2 transition (or moment) is described by the following one-body operator of nucle- 
ons, 

r(E2) = ^earfy(2)(f,)], (2) 

i 

where i is the index of nucleons constituting the nucleus. Within the framework of the 
present shell model, we are dealing only with the valence particles in the pf-shell. Hence 
the summation in Eq.(||) will be restricted to the nucleons outside the '^''Ca core, and the 
effective charge should be employed in order to incorporate the core polarization effect. 

In many cases the effective charges are fitted phenomenologically to experimental data 
under the assumption of the constancy for many transitions in one or several nuclei. In this 
paper, however, we derive the effective charges from a microscopic standpoint. In order 
to take into account the core polarization effect, we use the method of Sagawa-BrownP^ , 
which was shown to work well at least for collective transitions between lowest-lying states 
in the vicinity of doubly magic core. Then Cj in Eq.(0) is replaced by e'^J^{j,j'), where p 
is the subscript for distinguishing between protons and neutrons, and j {j') denotes the 
initial (final) single-particle orbit of the i-th nucleon. 

Although the Kuo-Brown interaction is calculated by assuming the harmonic oscillator 
single-particle wavef unctions, here we adopt the radial part of the single-particle wavefunc- 
tions given by the Hartree-Fock (HP) calculation. Although we abandon the consistency 
between them, this does not give rise to notable differences as shown later. The HP cal- 
culation is carried out with the SGII Skyrme interaction!^ for the ^^Ni core, assuming 
the full occupancy of the 0f7/2 orbit. The isoscalar (IS) and isovector (IV) quadrupole 
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Figure 23: 

Isoscalar and isovector quadrupole response function from the ^^Ni closed core, obtained from 
the HF+RPA calculation with the SGII Skyrme interaction. 

response functions are calculated by the random phase approximation (RPA)|]4T||. We 
searched giant quadrupole resonances (GQR) in the range of 10 < Ex < 70MeV, and 
find a single isolated IS-GQR peak at Ex ^ 17MeV and twenty-four IV-GQR peaks 
distributed over Ex = 20 ~ 35MeV. The response function is defined by 

S{E) = Y.\{u;i'^\T('^\0)mE - E^), (3) 

n 

where the state |0) is the HP state, Icu^'^^) an excited RPA state with the angular momen- 
tum A. In this case A = 2, and T^^^^ stands for the operator r'^Y^'^^ or r 

We show 

these response functions in Pig.^ taking into account the escaping width. The spreading 
width is not included in this calculation. By classifying the transition densities at the 
peak energies through their shapes, we select nine most prominent IV-GQR peaks as well 
as the IS-GQR peak. The transition strengths are assumed to concentrate in the selected 
GQR peaks. The adopted peaks and the corresponding transition strengths are listed in 
Table |. 

The core polarization effect caused by these GQR states is incorporated into the single- 
particle wavefunctions with the mixing amplitudes evaluated by the perturbation. When 
a single-particle state is denoted by the renormalized single-particle state is written 
as 

b) = |j) + E«n,'(j)k'^x/;j). (4) 
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Table 4: 

Excitation energies (MeV) and transition strengths (e^fm'^) of the GQR states adopted for the 
calculation of the core polarization effect. The transition strengths are summed over the states 
which possess similar shapes of transition densities. 



mode 


Ex 


BCE2) 


IS 


16.69 


800.9 


IV 


21.83 


29.3 




24.64 


93.4 




25.60 


38.7 




27.19 


9.4 




27.69 


13.2 




29.77 


35.2 




30.21 


31.3 




30.93 


67.5 




31.86 


284.3 



Here anj'(j) represents the mixing amplitude evaluated by the perturbation, 

(a;Wx/;j|V;.es|j) 
«nj'U) = : m— > (5) 

where Vres is the residual particle-hole interaction. Then the renormalized transition 
density is calculated by the perturbation as 

{]'\\T^'\r)\\j) c {f\\T^'\r)\\j) + Y: {anAjmT^'\r)\\u;W x f; j) 

n s, 

+a„,,(/)(a;Wxi;/||rW(r)||j)} 
= {f\\T^'\rm + X: |a„,.(j)^^(0||rW(r)||4^)) 

+(-)^+^-^'a„,(/)^|^(a;W||TW(r)||0)| . (6) 

In the RPA calculation and the evaluation of the mixing amplitudes, the residual interac- 
tion Vres is derived from the SGII Skyrme interaction, consistently with the HF calculation. 
After taking expectation values of momentum-dependent terms with respect to the HF 
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state, the Skyrme interaction Hsk only depends on the isoscalar nucleon density pis(r) 
and the isovector nucleon density pi\/{r). Therefore the residual interaction is written as 



1 S'Hsi, 



2 6p^,s 



1 S^H. 



Sk 



I t 



(7) 



2 Spjy 

where |o means the evaluation by the HF state. All steps of this renormalization are carried 
out at the ^^Ni core. It should be stressed that this calculation includes no additional 
adjustable parameters besides those contained in the Skyrme interaction. 

The reason why we start from the ^^Ni core is that the nuclei under consideration 
have masses closer to ^^Ni than to ^''Ca, as a whole. The difference in mass number 
affects (r^)s.p. (expectation values of for the single-particle orbits). It is noted that, 
from the ^^Ni core, we have quadrupole transition strength in the Ohu space. This Ohu 
strength should be handled principally within the shell model, since we have included the 
A; > configuration. The quadrupole transition strength within the Ohu space, however, 
seems to be exhausted in Ex < lOMeV, since hu ~ lOMeV at A = 56. Indeed we 
find only a single IS-GQR peak at Ex ~ 17MeV, which should be dominated by the 
2fiuj configuration. We search the RPA peaks for Ex > lOMeV, omitting the strength 
within the Ohu space. By this procedure, the E2 transition strength contributed by the 
k > 2 configurations, which stays within the Ohuj space but is neglected in the shell model 
calculation, cannot be included in the following calculation. It should be remembered 
that the E2 properties, particularly the collective ones, will provide us with a suitable 
test for the convergence with respect to the k <2 truncation. 

This renormalization procedure naturally introduces a j-dependence into the effective 
charges. These j-dependent effective charges are shown in Table ^. It is seen that the 
j-dependent effective charges are actually insensitive to j. This fact justifies the usual j- 
independent effective charges for the E2 transition. Recall that the j-independent charges 
were adopted in Ref.p], by adjusting the E2 transition probabilities in ^^Fe; e'^ = 1.4e 
and e*^ = 0.9e. Comparing the microscopic values in Table |^ with those j-independent 
effective charges, we have a good agreement for proton charges and a slight suppression for 
neutron charges. This implies that the core polarization effect is slightly underestimated 
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Table 5: 

E2 effective charges evaluated by the method of Sagawa-Brown. 



j 


f 




ef{j,j') 


Of'7 /O 

"-'^7/2 


"J- 7/2 


1.354 


0.677 


0f7/2 


Oh/2 


1.482 


0.841 


0f7/2 


lP3/2 


1.488 


0.788 


Oh/2 


Oh/2 


1.336 


0.643 


Oh/2 


lP3/2 


1.436 


0.712 


Oh/2 


lPl/2 


1.411 


0.690 


lP3/2 


lP3/2 


1.341 


0.601 


lP3/2 


lPl/2 


1.343 


0.604 



in both IS and IV modes. The influence of the k > 2 configurations within the pf-shell 
might be present in ^^Fe and this influence may account for the difference in the neutron 
effective charge. Note that the single-particle matrix elements of obtained by the HF 
calculation is very similar to those in the harmonic oscillator approximation, when we 
use the usual oscillator length b = yl^/®[fm] at A = 56. In practice, the ratios of the 
HF matrix elements to the harmonic oscillator ones are 0.93~1.05. Although there ought 
to be a certain decrease of the size of the single-particle wavefunctions as A becomes 
smaller, we will neglect this effect in the following calculation. This decrease amounts 
only to a 7% reduction of B(E2) at A — 50, by estimating it under the assumption of 
(r2),.p. oc 62 oc 
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4.2 Z = even, N = 28 nuclei 

In Table ^, the calculated E2 properties are compared with the measured ones, for Z = 
even, N = 28 nuclei. 

We first see B{E2] 2^ —>■ Of) in these nuclei. In the description only with the A; = 
configuration, the identical transition rates are predicted between ^°Ti and ^^Fe due to 
the particle-hole conjugation, as far as we neglect the change of the size of the 0fi/2 
single-particle wavefunction. The present k < 2 calculation inherits this feature to a 
certain extent. Provided that the seniority is conserved as discussed in the work of Horie- 
Ogawa^l, we also have the same B{E2) value for ^^Cr as for ^'^Ti or ^^Fe, within the k = 
space. In the present calculation, B{E2] 2^ —>■ Of) is distinctly enhanced in ^^Cr, due to 
seniority mixing]^] and the influence of the A; > configurations. 

On the other hand, the B{E2) values of ^^Fe is about twice as large as that of ^°Ti. 
The B{E2) value in ^^Cr is enhanced from that in ^°Ti, while it is close to the value 
in ^^Fe. Although the present calculation is far from reproducing the data for ^'^Ti and 
^^Fe, which has been a problem for a few decadesH, the result in ^^Cr yields a good 



agreement. As has been discussed in Subsection ^71|, the k < 2 truncation might be 
insufficient in reproducing some properties of ^^Fe. The k > 2 influence may account for 
the underestimate. In ^°Ti, the sd-shell contribution might be considerable. It is not easy 
to say whether the stronger mixing of the sd-shell configurations leads to an overestimate 
or an underestimate, because the adopted effective charges contain the excitation from the 
sd-shell. The success in ^^Cr, however, may indicate that the amount of core polarization 
is appropriately evaluated. Since effects of the k > 2 configurations are not included in 
calculating the effective charges, the wavefunctions look almost convergent by the k < 2 
truncation in this nucleus. It would be impossible to reproduce the rapid change of the 
B{E2) values without either more precise wavefunctions including the k > 2 configurations 
and the sd-shell excitations explicitly, or a strong particle number dependence of effective 
charges. 

In ^^Cr, there are two 4"'" states around 2.5MeV. Judging from i?(E2; 4"'" 2f) and 
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B{E2; Qf — > 4"^), the lower 4+ level in the calculation corresponds rather to the observed 
42 state, whereas the higher one to the observed 4^ state. This inversion has been 
taken into account in Fig.^ and Table |^. Through this assignment, we recognize a good 
agreement between the calculated and measured B{E2) values. 

In the k = space, the calculated quadrupole moments of the 2i states should have 
the same magnitudes in ^°Ti and ^^Fe, while the signs are opposite to each other: positive 
in ^°Ti and negative in ^^Fe. Although the signs are maintained, the magnitudes are 
different in the k < 2 calculation. The calculated quadrupole moment in ^^Fe is almost 
twice larger than that in ^°Ti. The experimental data contain too large errors to confirm 
this tendency. In ^^Cr, Q{2f) should be zero under the seniority conservation in the k = 
space. This obviously contradicts with the experiment, which gives a definitely negative 
value. The present calculation reproduces this property. 

The B(E2; 2i 0^) value of ^^Ni is already shown in Table ^ The result of the k < 4 
calculation, as well as that of the k < 2 calculation, are compared with the experimental 



value exhibited in Ref . ||33|| . Both the k <2 and k < A calculations give similar values, in 
good agreement with the datum. 
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Table 6: 

B(E2) values (e^fm^) or E2 static moments (efm^) in Z = even, N = 28 nuclei. The expression 
i (/) for the second (third) column denotes initial (final) state. The 'Cal.' values are obtained 



nucl. 


i 


/ 


Cal. 


Exp. 


50Ti 






88.0 


58. 


± 


9. 








86.8 


60. 


± 


13. 








40.1 


34.4 ± 


1.4 








1.6 


18. 


± 


9. 








9.0* 


8. 


± 


16. * 


52Cr 






102.7 


131. 




6. 








87.9 


84. 


± 


22. 








19.6 


53. 


± 


14. 




2+ 




0.05 


0.06± 


0.05 




2+ 

^2 


2t 


113.4 


150. 




35. 



nucl. 


i 


f 


Cal. 


Exp. 


52cr 






74.4 


59. ± 2. 




6^ 




13.2 


29. 7± 1.2 








0.5 


127. ± 8. 




3^ 




0.2 


7. ± 5. 




2^ 




-5.8* 


-14. ± 8. * 








73.3 


128.5± 4.8 




< 




68.6 






2t 


2t 


-16.4* 


-5. ±14. * 



*) Quadrupole moment. 
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4.3 Z = odd, N = 28 nuclei 

Table |^ shows the E2 quantities in Z = odd, N = 28 nuclei. 

As in the case of ^''Ti and ^^Fe, V and ^'^Mn are closely related through the particle- 
hole conjugation within the A; = space. If the seniority is a good quantum number, any 
corresponding transitions between low- lying states have the same B{E2) value in these 
two nuclei, while the quadrupole moments of the corresponding states have identical 
magnitudes but opposite signs. The present calculation maintains this nature for many 
E2 transitions. The transition from to is a typical example, which is also 
confirmed by the experiment. On the other hand, the most significant deviation is seen 
in B{E2; [|]^ [flT)- The i?(E2) value of this transition in ^-"^V is almost five times 
larger than that in ^^Mn. These properties are reproduced quite well by the present k < 2 
calculation. 

The results obtained in the present work are remarkably good in V. There are certain 
discrepancies in ^^Mn, which cannot be removed without discarding the particle-hole 
symmetry. This could be an infiuence of the k > 2 configurations. 
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Table 7: 



i?(E2) values (e^fm'^) or E2 static moments (efm^) in Z = odd, N = 28 nuclei. The 
experimental data are taken from Refs.p6|, p^, |30|, ^ |. 



nucl. 


i 




Cal. 


Exp. 


49Sc 


[|]r 


[|]r 


-20.9* 




Sly 


[fir 


[iir 


188.3 


169. ± 34. 




[fir 


[|]r 


75.6 


89. ± 10. 




[fir 


[fir 


101.1 


118. ± 16. 




mr 


[|]r 


95.2 


125. ± 21. 




[fir 


[iir 


31.6 


37. ± 6. 




[fir 
[f]r 
[|]r 


[fir 

mi 


29.9 
71.6 

-5.7* 


29. ± 5. 
66. ± 8. 
-4.3± 0.5* 



nucl. 


i 




Cal. 


Exp. 


53Mn 


[fir 


iBi 


180.5 


165.± 35. 




ir 


iBi 


62.0 


158.± 13. 




[fir 


[fir 


22.1 


18.± 6. 




mi 


iBi 


83.4 


151. ± 21. 






[|]r 


35.4 


83.± 11. 




[fir 


[fir 


33.9 


44.± 7. 






[fir 


0.1 


180. ±140. 




mi 


mi 


47.7 


64.± 10. 




iBi 


[|]r 


7.5* 




55 Co 


[|]r 


[|]r 


20.6* 





Quadrupole moment. 
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4.4 Z = even^ N = 29 nuclei 

We show the E2 properties oi Z = even, N = 29 nuclei in Table ^. 

Generally speaking, large errors in experiments make any precise comparison difficult. 
We only discuss the nearly degenerate | states. Let us consider B(E2; [|]]^ [flT) 
B{E2] ['^]2 [|]r)- The measured B(E2; [|]|f value points out the collectivity of 

this transition. In ^^Cr, the calculated transition strength from the higher | level agrees 
with the measured value much better than that from the lower | level. We therefore 
conclude that, in ^^Cr, the calculated higher | level should correspond to the observed 
state, whereas the lower one to the observed [|]2 state. This assignment has been 
taken into account in Fig.^and Table ^. The transition rates from [y]i and to these 
I states support this conclusion. The transition probabilities from the | states to [^]{ 
also prefer this reversed correspondence. On the other hand, in ^^Fe, as far as we can 
judge from the transitions to the ground states, the energy sequence of the calculated 
and [|]^ states seems to be correctly reproduced, as is assigned presently. 
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Table 8: 



i?(E2) values (e^fm'^) or E2 static moments (efm^) in Z = even, N = 29 nuclei. The 
experimental data are taken from Refs.p8|, |3^, ^ . 



nucl. 


i 




Cal. 


Exp. 


5lTi 




on 


102.5 










L2J1 


[Il7 
L2J1 


82.7 


930 


± 200. 






L2J1 


96.4 


3^0 


± 110. 






L2J1 


1.0 


49 


j- 


1 n 

iU. 






[^ir 


80.6 


1 00 


± 


21. 




L2J2 


[llr 


20.4 


382 


± 


79. 






L2J1 


37.5 


376 


± 


12. 




L 2 Jl 


[llr 


110.1 


95. 


± 


17. 




[fir 


[fir 


59.8 


62. 


± 


24. 




[fir 


9 1 


-11.3* 








53 


[2J1 




209.0 










L 2 J 1 


L2J1 


44.8 


169. 


± 


28. 






oh 

L 2 J 1 


146.1 


124. 


± 


11. 




[|]r 


1-2J1 


10.3 


76. 


± 


19. 




[|]2" 


1.2J1 


0.2 





26± 


0.06 




[|]2- 


[2-I1 


0.0003 


1 


9 ± 


1.7 




[i]2" 


9 1 


46.9 


43. 


± 


19. 




L ^ J ^ 


[|ir 


27.7 


470. 


± 470. 






[|ir 


132.9 


110. 


± 


12. 






12J2 


5.2 








nucl. 


i 


J 


Cal. 


Exp. 


53Cr 




[|ir 


18.2 












[II2- 


66.1 










ir 


[iir 


0.8 










ir 


[II2- 


306.8 


270 


±110. 




[i]2 


[fir 


21.9 


130 


±260. 




[fir 


[fir 


90.9 


30 


± 


12. 




[fir 


[fi2- 


356-1 


17 


± 


11. 




[fir 


[flr 


-15.4* 


-15 


± 


5. * 


err 











4.5 Z = odd, iV = 29 nuclei 

In Table ^ the E2 properties are shown for Z = odd, N = 29 nuclei. 

We can notice the overall agreement between the present calculation and the data, 
though disagreement in the order of magnitude is seen in B{E2; li 2^) in ^^Mn, 
S(E2;3^ ^ 3^), S(E2; 1^ ^ 3^) and S(E2;6^ ^ 5^) in ^^Co. As will be discussed 
later, the Ml properties of the relevant states are reproduced fairly well. Therefore the 
shell model wavefunctions of these states are probably adequate. These discrepancies may 
originate from the influence of the k > 2 configurations, or from the variation of the core 
polarization effect, which is possible since the initial states of the relevant transitions have 
rather high energy. 
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Table 9: 



i?(E2) values (e^fm^) or E2 static moments (efm^) in Z = odd, N = 29 nuclei. The 



experimental data are taken from Refs.|27, 31, 33 1. 



nucl. 


i 


/ 


Cal. 


Exp. 


50Sc 


"1 


"1 


24.0 






^1 




3.5 










-25.8* 




52y 


1 


3^ 


62.0 


< 9. X 10^ 




5+ 


3+ 


75.5 






1+ 
1 


3| 


96.3 






4^ 


3^ 


31.0 






3^ 


3^ 


3.8* 




54 Mn 


2^ 


3^ 


6.8 






4^ 


3^ 


83.9 


29. ± 11. 




5t 


3^ 


58.4 


109.± 




5t 


4| 


118.2 


280. ± 22a 




3t 


3^ 


0.4 


< 280. 




3+ 


2t 


172.6 


194.± 97. 




3+ 




90.8 


i20.± 1g: 








35.9 


5o.± ^iS: 








20.3 


17.± 




3+ 
"^3 


4^ 


31.6 


22. ± f?: 



nucl. 


i 


f 


Cal. 


Exp. 


54 Mn 


6^ 




0.003 


< 0.02 




1^ 


2^ 


2.4 


145. ±97. 




3^ 


3| 


31.6* 


33. ± 3. * 


56 Co 


3^ 


4^- 


47.1 


> 3.6 




5^ 


4^ 


58.3 


760. ± Z'o: 






4^ 


18.5 


< 70. 






3^ 


40.2 


< 52. 






4^ 


7.1 


25. 5± 8.9 




2^ 


3^ 


3.5 


9. ± 1 




5^ 


4^ 


8.7 


13. ± fo. 




3+ 

•^2 


4^ 


5.2 


10.2± 8.9 




3+ 
•^2 


3^ 


3.6 


89. ±51. 




3+ 




52.2 


382. ± 64. 






3| 


0.3 


25. ±13. 






5^ 


0.3 


< 32. 






5^ 


46.6 


0.5± li 






4^ 


22.3* 





Quadrupole moment. 
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4.6 Z = even, = 30 nuclei 



Table [T^ exhibits the i?(E2) values and quadrupole moments in Z = even, N = 30 nuclei. 

In, a good agreement is obtained. The only problem found in ^^Ti and ^^Cr is the 
transition from 2;^ to 2f in both nuclei. It should be noted that for the 7-transition from 

to 2i in ^^Cr, a different E2/M1 mixing ratio from the previous value has recently 



been reported [|T0[. The B(E2) value should be reduced if this new value is adopted. 



In ^^Fe, it is found that the present calculation tends to underestimate the B(E2) 
values. The k > 2 configurations will lead to expected larger collectivity. Keeping this in 
mind, we can recognize overall agreement also in this nucleus. The 6^ state has a large 
E2 transition probability to 4^, indicating a certain collectivity of the state. This B{E2) 
value is even larger than i?(E2; 6]*' Af). This feature is properly reproduced in the 
present calculation. There are certain discrepancies in B{E2; if — > 2^''), i?(E2;42' — > 4^^), 
i?(E2; Qf — » 4^) and -B(E2; 6^ — » 6f), whose initial states have rather high energy. 

The B{E2; 2^ Of) is quite underestimated in ^^Ni. It has been known that larger 
effective charges than those used in systematics are needed in order to adjust B{E2; 2f 
of) in ^^Ni0. Shimizu and Arima described the B{E2) value by taking into account the 
k = 1 and 2 configurations, based on the pseudo-SU(3) picture |Q. Although the present 
calculation includes the k = 1 and 2 configurations fully, anomalously larger effective 
charges are still necessary to reproduce the B{E2) value in ^^Ni, as in the case of ^"^Fe. 

As mentioned in Section ^, the calculated 0^ state of ^^Ni corresponds rather well in 
energy to the second lowest 0+ state observed so far. The B{E2) values to 2f and 2^ are 
consistent with this assignment. 
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Table 10: 

i?(E2) values (e^fm'^) or E2 static moments (efm^) in Z = even, N = 30 nuclei. The 



nucl. 


i 


f 

J 


Cal. 


Exp. 


52^- 


1 


1 


100.2 


1 38 + 


104. 
92. 




2+ 

2 


1 


13.7 


31 + 


23. 
14. 




2+ 


2| 

1 


53.7 








1 


2^ 
1 


134.0 








1 


4| 
1 


88.6 


123. lb 


22. 




i 


2| 


-7.5* 






54 Q J. 


1 


"1 


179.1 


173.6 ± 


3.0 




i 


2| 


237.6 


303 + 


97 




2+ 


1 


1.2 


10 9 + 


4 8 




2+ 


2| 
1 


71.9 


1 09 + 


36 




2+ 


2| 

1 


5.6 


~> 991 






"1 


1 


196.8 


Z lO. ZIZ 


fil 




2^ 
1 


1 


-24.5* 


—21. ± 


8. * 


56pg 


2^ 

1 


1 


164.8 


91 q Q _|_ 






4i^ 
1 


2| 
1 


229.6 


305 + 


64 




2+ 

2 


1 


2.7 


6.4 ± 


7.6 
6.4 




22^ 


2| 


4.5 


31. ± 


14. 




0+ 




6.7 


31. ± 


15. 




2+ 




0.6 


0.97± 


0.11 




2^ 




13.7 


16.5 ± 


5.1 






2t 


0.0004 


110.7 ± 


8.9 




42^ 


2| 


0.0008 


1.0 ± 


1.1 

1.0 




2+ 

^4 


0^ 


17.2 


10.2 ± 


5.1 




2| 




18.0 


15. ± 


19. 
15. 




6^ 




33.4 


50.9 ± 


5.1 



nucl. 


i 


/ 


Cal. 


Exp. 


56pg 




2^ 


13.0 


7.6 


± 6.4 




3^ 


4^ 




10.2 


± 3.8 




3^ 


2t 


149.3 


430. 


±140. 



4.7 Z = odd, N = 30 nuclei 

In Table 11, the calculated and measured E2 properties^] are compared for Z = odd, 
N = 30 nuclei. 

The agreement of the calculation with the data is good in ^^V and ^^Mn. Although 
the data in ^^Co are not necessarily reproduced well, we point out that the orders of 
magnitude are correct for most transitions. 



^ The (5([|]i ) value of ^^Co quoted in Ref. jsj] is questionable, since it is not even referred to in the 
original paper |43|. 
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Table 11: 



i?(E2) values (e^fm^) or E2 static moments (efm^) in Z = odd, = 30 nuclei. The 
ex perimental data are taken from Refs.pS], |3^, 3^ . 



nucl. 


i 


f 


Cal. 


Exp. 


^^Sc 


[ill 




50.0 






rill — 
[¥]l 


r 7i — 


22.7 






r9i- 


r7i- 


28.2 






r5i- 


r7i- 


3.4 






r7i- 


r7i- 


—23.8* 






r5n — 
[fll 


r7n — 


253.2 






[31- 
i2Jl 


[71- 
Wl 


ibz.z 


i4o.± iz. 




mi 




153.3 


183. ± 28. 




r 1 5 T — 

mi 


mi 


157.7 


154. ± 35. 




r7n — 
[|]l 




-10.7* 




55l\ yT 

iVin 


r7i- 
[2]! 




287.2 


335. ± 99. 




[fll 


[fll 


85.8 


109. ± 14. 




rQi — 
[fll 


r7i — 


195.6 


286. ± 37. 




rill — 

mi 


r7i — 


152.6 


175. ± 15. 




mi 




160.3 


< 87. 




[-2]! 




36.2 


144.± 12. 








19.6 


41. ± 26. 




[f]r 


ir 


148.5 


99.± 50. 




[f]r 




89.9 


40. ± 11. 




mi 




164.2 


96. ± 22. 




mi 


mi 


52.0 


32. ± 20. 




ir 




33.5* 


33. ± 1. * 


nucl. 


i 


/ 


Cal. 


Exp. 


s^Co 


[|]r 


[|]r 


135.3 


241. ± 29. 




[fir 


[|]r 


0.08 


6.0± 1.3 






[fir 


480.8 


< 23. 




mr 


[|]r 


34.1 


78.2± 7.8 






[fir . 


q114.9 


370. ±100. 






fCi-. 


4.9 


0.1± 1.6 



5 Ml properties 



5.1 Single-particle parameters 

The Ml properties are usually described by the following one-body operator, 



r(Mi) 

For the free nucleon, we have 



■ h + 9s,T^i e proton) 
/Xi = <( , (9) 

Ss^u^i V' ^ neutron) 

in units of /iN, where U (sj) denotes the orbital (spin) angular momentum of the i-th 
nucleon, and the are spin (^-factors with g^^'^^ = 5.58, g^Jl'^ = —3.82. It is known, 
however, that the core polarization and the meson exchange current contribute to the Ml 
properties. The former appears to be the price for the truncation of the model space in 
the shell model calculation. Especially the quenching of the spin (^-factor is required in 
most cases. Considering the orbit dependence of the single-particle parameters as in the 
E2 case, we can parametrize the Ml operator as 

A, = gf^iniyh + gt{nl)h + g^M^. n'l')[Y^'\U)s,f~^ (p = vr, z/), (10) 

where the gis are orbital ^'-factors. The (^-factors are, in principle, dependent on n and 
/ quantum numbers of the single-particle orbit. We adopt the microscopic parameters 
that Towner calculated with the single-particle wavefunctions in the harmonic oscillator 



approximation 1 44 ] 



As in the calculation of energy levels and E2 properties, the mass-number dependence 
of the single-particle parameters is neglected. In the harmonic oscillator approximation, 
only the radial part varies with mass number. The matrix elements of the / and s opera- 
tors, which form the bare Ml operator, do not depend on the radial part. Only a part of 
the correction coming from the core polarization and meson exchange effects will change 
with A. Therefore the v4-dependence of the Ml single-particle parameters will be weak. 
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Table 12: 

B(y[l) values (/i^) ^1 static moments ()Un) in Z = even, N = 28 nuclei. The 'Cal.' values 
are obtained by the present calculation. The experimental data (Exp.) are taken from 



Refs. 



nucl. 


i 


/ 


Cal. 


Exp. 


50xi 






2.388* 


2.4 ± 0.8 * 


52Cr 






0.001 


0.0006± 0.0001 








0.009 


0.082 ± 0.027 




4 




0.920 






4^ 


4 


0.216 


0.358 ± 0.072 



nucl. 


i 


/ 


Cal. 


Exp. 


52Cr 


3| 




0.009 


0.013±0.005 




2^ 




2.732* 


3.00 ±0.50* 


54pg 


2^ 




2.732* 


2.4 ±0.3 * 




6| 


6^ 


8.334* 


8.22 ±0.18* 



Magnetic moment. 



5.2 Z = even, N = 28 nuclei 



The calculated Ml properties are shown in Table |T2], in comparison with the experimental 
data, for Z = even, N = 28 nuclei. 

All the measured magnetic moments are reproduced very well in ^°Ti, ^^Cr and ^^Fe. 
In the k = configuration, the particle-hole conjugation leads to the same fi{2i) between 
^°Ti and ^^Fe. This situation approximately holds both in the experiment and in the k < 2 
calculation. The contributions of the higher configurations, which have been suggested 
with respect to the E2 quantities, are not clear in the magnetic moments. The fi{2f) 
value of ^^Cr is close to those of ^°Ti and ^^Fe, as is expected if the seniority conservation 
is satisfied in the k = space. It should be mentioned that a smaller experimental value 
of fi{2t) = 2.00 ± 0.32[^^] has been indicated for ^^Fe recentlygB . 



In ^^Cr, the B{M1) values from 4J^ imply that the Af and 4^ assignment performed 
in Section H is reasonable. 
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5.3 Z = odd. N = 2S nuclei 



The Ml quantities in Z = odd, N = 28 nuclei are displayed in Table [13. 

All the measured magnetic moments are reproduced remarkably well, highlightening 
the adequacy of the present description with respect to the proton degrees of freedom. 
According to the particle-hole symmetry in the k = space, the same magnetic moments 
are expected between ^^V and ^^Mn, as well as between ^^Sc and ^^Co. This equality is 
satisfied approximately. 

As in the case of ^^Cr, the qualitative trend of the B{M1) values is reproduced in ^^V 
and ^^Mn. From a quantitative point of view, we point out that the I? (Ml) values are 
generally underestimated in the Z = odd or even, N = 28 nuclei. 
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Table 13: 

-B(Ml) values (/x^) or Ml static moments (^n) in Z = odd, N = 28 nuclei. The experimental 



data are 


^aken from Refs. 


2| 




nucl. 


i 


/ 


Cal. 


Exp. 


49Sc 


Sr 


Sr 


5.152* 




Sly 


[|]r 




Sr 


0.0003 


0.0053 ± 0.0003 




[iir 




ir 


0.0008 


0.00007± 0.00002 




[|]r 




Sr 


0.0005 


0.0006 ± 0.0002 




[|]r 




¥]r 


0.033 


0.082 ± 0.045 




Sr 






5.008* 


5.149 ± 0.000* 




[fir 




[fir 


3.457* 


3.86 ± 0.33 * 


nucl. 


i 


/ 


Cal. 


Exp. 


53Mn 


[fir 




Sr 




0.0001 


0.0046±0.0004 




[fir 




[fir 




0.0009 


0.041 ±0.005 




[fir 




©r 




0.00006 


0.0021±0.0005 




[fir 




[fir 




0.0009 


0.0026±0.0003 




mi 




[fir 




0.024 




0.376 ±0.072 




iBi 




[|ir 




4.996* 




5.024 ±0.007* 




[fir 




[fir 




3.488* 




3.25 ±0.30 * 


55 Co 


[|]r 


[|ir 


4.966* 


4.822 ±0.003* 



*) Magnetic moment. 
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5.4 Z = even, = 29 nuclei 



Table shows the Ml properties for Z = even, N = 29 nuclei. 

We first look at ^^Ni, which could be a good test of the neutron degrees of freedom. 
In Table 0, the sign of the experimental magnetic moment of , which has not been 
specified|^, is conjectured from the calculated value. The adopted sign is consistent 



with the Schmidt value. The /i([|]r) ^-iid B(M1] [|]r) values are overestimated in 

the present calculation. Since the neutron orbital angular momentum hardly affects the 
Ml observables, it is suggested that further quenching of gs,u is required for more precise 
description. Though within the k = configuration the Ml transition from [|]]^ to 
[|]|f is possible only through the [y'^^^sJ^^-'-term, the measured probability is significantly 
larger than expected from a reasonable gp^u value. Analogously, the present calculation 
underestimates B(M1] by several orders of magnitude. In order to remedy 

this problem, we should include higher configuration into the wavefunctions, or two-body 
terms into the Ml operator. 

It should be noticed that the calculated /i([|]r) and /u([|]2 ) of ^^Cr have opposite signs 
to each other. The observed prefers to be positive. This fact is consistent with 

the present assignment based on the B{E2) values in Section ^, in which the indices of 
the calculated two | states have been inverted against their energy sequence. The signs 
of the Ml moments of the lowest two | states are correctly reproduced also in ^^Fe, even 
though the absolute value of /i([|]2^) is underestimated to a certain extent. 

The calculated Ml transition rates agree with the data, at least within the order of 
magnitude, for ^^Ti, ^^Cr and ^^Fe. The only exception is B(M1; [^]2 in ^^Fe, 

which is not a serious problem since this B{M1) value is small both theoretically and 
experiment ally. 
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Table 14: 

B{y[l) values (/U^) or Ml static moments (/^n) in Z = even, N = 29 nuclei. The experimental 



data are 


taken from Refs.p^ 


, 30,^,11. 


m 1 r*! 




f 
J 




Exp. 


i 1 


fir 


l2h 




0.734 ± 0.090 




[2J1 


[2J1 


0.049 


0.059 ± 0.018 




L2J2 


[2J1 


0.002 


0.016 ± 0.003 




[2J2 


[2J1 


0.068 


0.013 ± 0.003 




L2J2 


I2J1 


0.146 


0.170 ± 0.032 




L2-I2 


12J1 


0.085 


0.127 ± 0.004 




L2J2 


L2J1 


0.146 


0.516 ± 0.016 




L2J2 


12J1 


1.739 


3.22 ± 0.11 




[fir 


[fir 


-0.989* 




53 


I2J1 


ir 

I2J1 


0.653 
0.003 


0.093 ± 0.013 




L2J1 


L2J1 


0.273 


0.082 ± 0.016 




\^]~ 
L 2 


1 

L 2 -1 i 


0.001 


0.0075± 0.0005 




f^r 

L 2 J2 


L 2 J i 


0.00006 


0.029 ± 0.003 




f^r 

L 2 J2 


1 

L 2 -1 i 


0.076 


0.050 ± 0.016 






TT 1 

L2-I i 


0.310 


0.269 ± 0.090 




9 
oh 

L 2 J -L 

oh 

L 2 -1 -L 


9 

L 2 -1 J- 

f^r 

9 

L 2 -1 ^ 


0.00000 

0.001 

0.618 


0.315 ± 0.045 


nucl. 


i 


/ 


Cal. 


Exp. 


53cr 




[fir 


0.533 


3.6 ±7.2 




[fir 


[fir 


-0.442* 


-0.475 ± * 




[|]r 


[|]r 


2.435* 


2.8 ±4.9 * 




©2 


©2 


-0.719* 




55 pg 


[^]r 


[fir 


0.844 






[fir 


[fir 


0.0004 


0.005 ±0.002 




[|]r 


[fir 


0.123 


0.023 ±0.016 




[|]2 

[fir 


[fir 

r- '7-, 


Q^.^0007 
0.317 


0.0050±0.0004 
0.093 ±0.027 



5.5 Z = odd. N = 29 nuclei 



Table |T5| shows the Ml properties for Z = odd, N = 29 nuclei. We obtain good agreement, 
except for the extremely large /i(4^) and fJ-i^f) of ^^Mn which are beyond the ordinary 
shell model description. 
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Table 15: 



B{M.l) values (/x^) static moments (^n) in Z = odd, N = 29 nuclei. The experimental 
data are taken from Refs.(2^, 1^, ||, H]. 



nucl. 


i 


f 


Cal. 


Exp. 


nucl. 


i 


f 


Cal. 


Exp. 


50Sc 


"1 


1 


2.099 






1+ 


2t 
1 


0.030 


0.025± 0.014 




1 


"1 


1.748 






3^ 
1 


3| 
1 


2.987* 


3.282± 0.001* 




i 


1 


3.980* 






4| 
1 


4| 
1 


3.749* 


7.3 ± 1.4 * 


52v 




3^ 


1.528 


1.33 ± 0.50 




"1 


"1 


4.282* 


55. ±30. * 




1+ 
1 


1 


0.226 






3+ 

2 


3+ 

2 


4.226* 


< 18 




1 


"1 


0.110 




^^^Co 


3^ 


4^ 
1 


2.190 


> 0.10 






5t 


0.572 








4^ 


0.552 


0.72 ± 0.18 




3^ 


3| 


3.049* 






4^ 


4^ 


0.003 


< 0.010 


54Mn 


2t 


3^ 


1.158 








3^ 


0.035 


< 0.057 




4^ 


3^ 


0.100 


0.055 ± 0.005 






5| 


0.003 


< 0.021 




5^ 


4^ 


0.538 


0.591 ± 0.072 






3^ 


0.855 


0.97 ± 0.27 




3+ 


3| 


0.099 


0.116 ± 0.025 




5+ 

^2 


4^ 


0.223 


0.088± 0.036 








0.034 


0.052 ± 0.013 




3+ 

-32 


< 


0.103 


0.125± 0.072 




3+ 




1.182 


1.00 ± 0.22 




3+ 

"^2 


3^ 


0.006 


0.006± 0.003 




4 




0.499 


0.322 ± 0.072 




3+ 




0.727 


1.07 ± 0.54 




3+ 




0.420 


0.59 ± 0.14 








0.00006 


0.143± 0.072 




3^ 


4| 


0.245 


0.63 ± 0.14 




6^ 




0.142 


0.125± 0.072 




6^ 


5^ 


0.0004 


0.0005± 0.0002 




4^ 


4 


3.654* 


3.830± 0.015* 



*) Magnetic moment. 
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5.6 Z = even, = 30 nuclei 



The calculated Ml quantities are compared with the measured ones in Table 0, for 
Z = even, = 30 nuclei. 

Although the agreement of the B{M1) values is not necessarily good in ^^Fe, the 
orders of magnitudes are correctly reproduced for most transitions. As will be discussed 
below, the further quenching of gs,p improves the results considerably. The observed 
B{M1; if of) value is quite small. This nature is reproduced well in the present 
calculation. The discrepancy in B{Ml;2f — ^ 2^*^) is not improved by the Qs quenching. 
The cancellation among the elements of the density matrix brings this small B{M1) value. 
As discussed in Ref. [Q, the 23 state is a fully non-collective state. We should note that 
non-collective transitions are generally sensitive to details of wavefunctions. The {p,p') 
result confirms the adequacy of the present shell model wavefunction of this state, as 
will be mentioned in Section 0. 

In Ref.0], the n and / dependences in Eq.(|lOD were neglected, as is often assumed. 
Furthermore, Qp^^ = was assumed. Namely, the employed Ml operator was 

TiMl) = ^/^ E (91% + C^p)^ (11) 

where 

Lp = Y,h, Sp = J2s,. (12) 

The single-particle g-iactois of gf^ = 1.0, g'(^ = 0.0, g^^^ = O.bg^^^ and g^^^^ = 0.5g^Jl^ were 
adopted in Ref.P] so as to fit the data in ^^Fe, though microscopic calculations predict the 
quenching factor for g^p to be much closer to the unity |]4^, It was shown in Ref.|^ 
that by this phenomenological set of parameters the i?(Ml) valued in ^^Fe are reproduced 
reasonably well. 

It is found that the predicted 13" state of ^^Fe, as well as the 1^ state, has a relatively 
large B{M1) value to the ground state. This fact suggests that these two states contain a 
considerable fraction of a mixed-symmetry l"*" component. A recent experiment has also 
reported a mixed-symmetry 1^ strength around Ex ~ 3.5MeV^^. We will discuss this 
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Table 16: 

B{y[l) values (/U^) or Ml static moments (/^n) in Z = even, N = 30 nuclei. The experimental 



nucl. 


i 


/ 


Cal. 


Exp. 


52rpi 


2+ 


2^ 


0.456 


0.56 ± 






2^ 


0.610 


> 0.16 






2^ 


0.666* 


— 


54 Cr 


2+ 


2^ 


0.021 


0.023± 0.009 






2^ 


0.519 


> 0.10 




2^ 




1.467* 


1.12 ± 0.20* 


56 pg 


2+ 

^2 


2^ 


0.476 


0.233± 0.072 




2^ 


2| 


0.0005 


0.070± 0.009 




1^ 


0^ 


0.0004 


< 0.00001 




1^ 


2^ 


0.090 


0.040± 0.003 




4^ 


4^ 


0.841 


0.206± 0.038 




2^ 


2^ 


0.202 


0.109± 0.045 




31 


2^ 


0.028 


0.054± 0.013 




3^- 


4^ 


0.030 


0.098± 0.023 




3^ 


2+ 


0.026 


0.018± 0.007 



XL Refs. 


1, 31, P, 






nucl. 


i 


f 


Cal. 


Exp. 


56pg 


1+ 

^2 


0+ 


0.061 


0.047 ±0.013 




it 
^2 


2+ 


0.048 


0.054 ±0.014 




it 
^3 


of 


0.151 






it 
^3 


2t 


0.155 






2+ 


2+ 


0.012 


0.003 ±0.002 




6t 


6t 


0.274 


1.09 ±0.18 




2+ 


2+ 


1.379* 


1.20 ±0.20* 


58Ni 


2+ 

^2 


2+ 


0.156 


0.018 ±0.009 








0.00004 


0.0014±0.0003 




0+ 




0.0001 






03^ 




0.0005 


0.158 ±0.009 




2+ 
^3 


2^ 


0.018 


0.134 ±0.021 




2+ 
^3 




0.140 


0.555 ±0.090 




2| 




0.030 


0.070 ±0.029 




3^ 




0.308 


0.16 ±0.14 




2^ 


2t 


-0.168* 





Magnetic moment. 



point in more detail in a forthcoming paper. These two states are almost degenerate in 
energy {Ex ~ 3.5MeV) according to the calculation, while only a single state has been 
experimentally confirmed ||33||. It is desired to search for the other state. 

We do not have a good agreement for ^^Ni. This disagreement might imply a non- 
negligible infiuence of the k > 2 configurations. 
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5.7 Z = odd. iV = 30 nuclei 



The Ml properties are shown in Table [T^, for Z = odd, N = 30 nuclei. 

The Ml quantities in ^^Mn are reproduced remarkably well. Together with the agree- 
ment in the E2 properties, this fact suggests good convergence of the wavefunctions in 
this nucleus. The agreement in ^^Co is only for the orders of magnitude. 
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Table 17: 



B(Ml) values (/x^) or Ml static moments (^n) in Z = odd, N = 30 nuclei. The experimental 

data are taken from Refs.|2|, g^, Igg]. 



nucl. 


i 


f 


Cal. 


Exp. 


s^Sc 


L2 -1 1 


L2 J 1 


5.061* 




53v 


L2 J 1 


L2 J 1 


0.065 


> 0.027 




L2J1 


L2 J 1 


0.00001 


0.0032± 0.0003 




L2J1 


L2J1 


4.670* 




55Mn 




L2 J 1 


0.042 


0.076 ± 0.003 




L2 J 1 


L2 J 1 


0.085 


0.197 ± 0.021 






[|]r 


0.227 


0.224 ± 0.023 




[fir 


[fir 


0.147 


0.167 ± 0.036 




[fir 


mi 


0.154 


0.233 ± 0.054 




mi 


mi 


0.257 


0.251 ± 0.072 




[fir 


[fir 


3.378* 


3.453 ± 0.001* 








4.477* 


4.4 ±0.7* 


s^Co 


ir 


Sr 


0.304 


0.372 ± 0.036 






[fir 


0.652 


0.090 ± 0.011 




mi 


[|]r 


0.722 


0.877 ± 0.090 



nucl. 


i 




Cal. 


Exp. 


57Co 


[|]2- 


[|]r 


0.002 


0.022±0.003 






[fir 


0.115 


0.662±0.090 




[fir 


[|]r 


0.044 


0.236±0.032 




[l]2 


[|]r 


0.002 






iBi 


[|]r 


0.032 


0.003±0.001 




[|]3 


[fir 


0.775 


0.100±0.020 




[i]2- 


©r 


0.020 


< 0.034 






[fir 


0.028 


< 0.054 




[i]2- 


[fir 


0.001 


< 0.21 






[|]2- 


0.004 


0.36 ±0.13 




mi 


[fir 


0.718 


0.73 ±0.18 




mi 


[fir 


.#•730 


0.233±0.090 




[|]r 


[|]r 


4.729* 


4.727±0.009* 



5.8 Summary of electromagnetic properties 

We summarize the consequences for the electromagnetic properties. The E2 transition 
probabihties and moments are comprehensively reproduced, together with the effective 
charges derived by the Sagawa-Brown method. The good agreement with the measure- 
ment in the collective {i.e., enhanced) E2 transitions suggests that the convergence of the 
wavefunctions is reasonable in the framework of the k <2 truncation. We also acquire an 
overall reproduction of the Ml transition rates and moments, by using Towner's single- 
particle g'-factors. Particularly, most of the measured magnetic moments are precisely 
reproduced. The remaining discrepancies may be ascribed to the k > 2 configuration or 
the sd-shell contribution. 



6 Electron scattering form factors 

We calculate longitudinal form factors of inelastic electron scattering with an angular 
momentum transfer of two units (C2), from the ground state to several 2"*" states for 
a few even-even nuclei. The form factors provide us with richer information than the 
7-transition probabilities, owing to the off-shell photon exchange. We here attempt to 
see how well the present calculation reproduces the quadrupole collective feature. The 
method of Sagawa-Brown ppj] , which has already been sketched in Section ^ is used again. 
This method was shown to work very well for longitudinal form factors with collectivity 
for nuclei around the doubly magic core. 

The C2 form factors are calculated from the shell model density matrices, together with 
the renormalized single-particle transition densities of Eq. . The nucleon finite-size effect 



is incorporated in the dipole approximation! 49], and the contribution of the center-of-mass 



motion is subtracted in the harmonic oscillator approximation ||50| . The plane- wave Born 
approximation (PWBA) is employed, taking into account the Coulomb distortion effect 
in terms of the effective momentum transfer ges- It should be emphasized that this 
calculation includes no additional adjustable parameters, as pointed out in Section ||. 
The results for ^^Fe were already shown in Ref.P]. Here we mention that the results 
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Figure 24: 

(e, e') form factors from the ground state to the lowest four 2^ states of ^"^Cr. The sohd Hues 
show the sheU model results with the core polarization effect evaluated by the HF+RPA. Those 
without core polarization are shown by dotted lines. The crosses exhibit the experimental data 
taken from Ref . p2| . The transverse mode is not separated in the experiment. 



of this calculation agree with the data||51| very well. Particularly, the form factor in the 
excitation to 2^ is in a precise agreement. It is also remarkable that the collectivity of 
the 24 is correctly reproduced. 



The C2 form factors in Cr are shown in Fig.|2^. In the experiments the transverse 



mode was not separated, though most of the reaction is known to be dominated by the 



C2 mode[Q. The experimental data are described precisely with respect to the 2^^ state. 
As for the 22 state, the shape is reproduced well, though the absolute value is overall 
underestimated by a factor of 2. It is pointed out that this transition is not collective. 
There are no reported data with respect to the 2'^ and 2^ states. We show them, however, 
since they are good candidates for the mixed-symmetry state [pj . 

The form factors from to 2^ in ^^Fe and ^^Ni are shown in Fig.|2^ and Fig.p^, 
respectively. The absolute value is underestimated, and this is consistent with the B{E2) 
values discussed in Section This might indicate a worse convergence of the wavefunc- 
tions for these nuclei. Nevertheless, the shape of the form factors, especially the position 
of the peaks and the dips are reproduced very well. 
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Figure 25: 

(e, e') form factor from the ground state to the 2^ state of ^^Fe. See the caption of Fig. 24. The 

experimental data are taken from Ref . [^4|| . 



Figure 26: 



(e, e') form factor from the ground state to the 2f state of ^^Ni. See the caption of Fig. 24. The 

experimental data are taken from Ref. [p5| . 

7 Proton scattering cross sections 



Recently {p,p') experiments are performed for ^^Fe and ^^Cr|]T2|, with an incident proton 
energy of 65MeV. The differential cross sections and the analyzing powers are extracted. 
The data are analyzed in the distorted- wave Born approximation (DWBA), by using the 
present shell model density matrices. This DWBA calculation is carried out by Takamatsu 
and his collaborators [1^. We briefly review their results in this section. 



The Bonn-Jiilich effective interact ion |^6[ is employed as the interaction between the 
incident or scattered proton and the nucleons in the target nucleus. Although the M3Y 
interaction [57 1 and the Paris-Hamburg interaction! 5^ are tried also, the results were es- 
sentially the same. In order to reproduce the experimental data, an overall normalization 
factor is used for each transition. This factor would correspond to the core polarization 
effect. In fact, the adjusted normalization factors are consistent with the ones expected 
from the effective charges employed in the description of the B{E2) values. 
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The absolute value of the cross section in the excitation to the 2^ state is smaller by 
an order of magnitude than to the 22 or the 2f states in the experimental data of ^^Fe. 
This is consistent with the present calculation. In addition, the excitation to the 2^ state 
shows a strikingly anomalous angular distribution. The present shell model wavefunctions 
reproduce this anomaly remarkably well. 

The experimental differential cross sections are also excellently reproduced for ^^Cr|l^ . 



We have not seen any notable contradiction between experiments and the present 
calculation with respect to (e, e') and {p,p') data. 

8 Summary 

The properties of the low-lying states in 20 < Z < 28, 28 < < 30 nuclei have been 
investigated from a microscopic standpoint. The model space has been restricted to the 
pf-shell, while we have incorporated the excitation from 0f7/2 to (0f5/2lp3/2lPi/2) 
two particles. This space leads to one of the largest-scale shell model calculations that 
have ever been carried out. 

We have adopted the Kuo-Brown interaction on top of the '^'^Ca core, including the 
3p-lh correction. This effective interaction is derived from the Hamada- Johnston A^A^ 
potential, through the G-matrix calculation. Despite the failure near the beginning of the 
pf-shell, it has been shown for the first time, by this work, that this microscopic interaction 
produces an excellent description of the nuclei around the middle of the pf-shell. 

The observed energy levels are reproduced by the present shell model calculation in 
a wide range of energy, for any of the 20 < Z < 28, 28 < < 30 nuclei. The energy 
range of good agreement is typically Ex < 4MeV for even-even nuclei. Ex < 2.5MeV 
for odd-mass nuclei and Ex <, 2MeV for odd-odd nuclei, with the discrepancies of 6E <^ 
0.3MeV. These energy ranges are notably wider than those obtained in previous shell 
model studies[^, ^ 0, covering the region where candidates for the mixed-symmetry 



states were reported. Among a number of states dominated by the k = configuration, 
more than twenty states in the low energy region are found to have k = 1 probabilities 
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larger than k — ones. The k — 1 dominance of those states is consistent with the fact 
that the Horie-Ogawa calculation is not capable of describing them. Those states are also 
reproduced as well as the A; = dominant ones. The present calculation is undoubtedly 
more successful than any other calculation that has ever been reported in this mass region. 

In order to study the E2 properties, we have calculated the effective charges by ap- 
plying Sagawa and Brown's HF+RPA method. The measured E2 moments and transi- 
tion rates are reproduced fairly well, by the shell model wavefunctions together with the 
microscopic effective charges. The Ml properties have also been investigated by using 
Towner's microscopic parameters. This approach gives an overall agreement with exper- 
iment. There remains, however, some discrepancies in the electromagnetic properties. 
Most of the discrepancies seem to be due to the influence of the A; > 2 or the sd-shell 
configurations. The overall reproduction of the electromagnetic properties has conflrmed 
the reliability of the shell model wavefunctions. 

The (e, e') C2 form factors have also been calculated in the same framework as the 
E2 properties. Together with the {p,p') results, we have confirmed that the quadrupole 
collective features are successfully described by the present shell model wavefunctions. 

We have investigated the low-lying states of the nuclei in the middle of the pf-shell, 
from a microscopic standpoint. The present effective interaction is derived on a realistic 
basis, without fitting two-body matrix elements to observed energy levels, although there 
remains an uncertainty in the G-matrix calculation. By applying this interaction, some 
of real situations of the low-lying states of the middle pf-shell nuclei are successfully 
reproduced. It is remarked that the large-scale calculation including the excitation from 
0f7/2 to (0f5/2lp3/2lPi/2) ^wo uuclcons is crucial for this agreement. The successful 

results of this study enable us to place confidence, to a certain extent, in the ways of 
understanding the nuclear many-body system from the nucleonic degrees of freedom. 

The authors arc grateful to Prof. H. Sagawa for valuable discussion of the HF+RPA 
calculation. The authors thank Prof. A. Gelberg for careful reading the manuscript. The 
computers used for the numerical calculations are HITAC-S820/80 in Computer Center, 
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University of Tokyo, HITAC-S820/80 in Kanagawa Factory, Hitachi Corporation, and 
VAX6440 in Meson Science Laboratory, University of Tokyo. This work is financially 
supported, in part, by Research Center for Nuclear Physics, Osaka University. 
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